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type of reinforcement, constituent material properties, and the physical location and orientation 
of the reinforcing members. In the case of a textile reinforced composite, quantification of the 
geometrical properties of the reinforcement is rather complex. Thus, the first step in a detailed 
modeling approach is the development of methodologies for predicting and reporting the positions 
and orientation of the yarns contained within the RYE. 
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determined by fabric design. 



The yarns forming a unit cell are considered to be elastic bodies interacting with one another 
and subject to external loads. The center line of every yarn is represented as a Bezier curve 
interpolating a set of discrete support points. In this study, a continuous Bezier with degree n — 1 
(over n support points per yarn) was used. The location of the support points fully determines the 
current mechanical state of the unit cell. The yarns are assumed to have elliptical cross sectional 
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The Lagrangian (1) does not account for the dissipation of energy in the system. Of course, 
in the general case friction can be one of the major factors influencing the mechanics of the unit 
cell. It is important to note that there are no principal difficulties in introducing the corresponding 
friction term into the model. 



The center lines of yarns are represented in the model by Bezier curves [7] interpolating sets 
of support points. In addition, a number of intermediate points placed on the center lines of yarns 
are selected. At any of these intermediate points tangent, normal and bi- normal vectors are built 
[8]. Ellipses located in the planes perpendicular to the center lines and containing intermediate 
points as their centers represent the surfaces of yarns. 
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to a local extremum which can result in the incorrect prediction of elastic properties of the fabric? 
In order to increase the probability of finding the state of the system corresponding to the global 
minimum of L, a number of techniques are known, but none of them can guarantee the convergence 
to the global minimum. For the problem under consideration the following method proved to be 
effective: the positions of several support points of different yarns are varied simultaneously during 



a single step. In particular, two contacting yarns can be displaced together without penetrating 



the specimen are described in Table 1. 


The parameters of the fabric (dimensions of yarns, distance between the yarns) and the 
results of testing, characterizing elastic properties of warp and weft yarns (C*), have been taken as 




the energy of tension for all yarns) the elongation of specimen in the warp direction occurs mostly 
due to the crimp interchange between warp and weft yarns, the length of warp yarns being almost 
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During the minimization of L the out-of-plane (z) coordinates of the end points of both 
weft yarns and both warp yarns were kept unchanged. This allowed simulation of real boundary 
conditions of the experiment. The Poisson’s effect in the specimen, neglected in this formulation of 
boundary conditions, is believed to be small because of the small distance between the jaws of the 


Both experimental and analytical load-strain curves are shown in Figure 3 The values of 
applied external forces have been normalized by a number of warp yarns in the specimen. Good 
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Table 1: Para meters of Plain Weave Fabric Unde r Consideration 
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Figure 1: Graphical Rendering of Plain Weave Fabric With Nominal External Load (25 g) at 
Equilibrium 
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Figure 3: Comparison of Experimental and Predicted Axial Stress-Strain Response of Plain Weave 



Densified Fabrics 



NUMERICAL GEOMETRIC RESULTS 


The analysis was performed for the four specific braids that are currently being tested as part of 
the overall NASA activity. These braids are characterized according to Table 3: 
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Figure 4: Geometry of Triaxial Braid Model. LSS, Front View. 



























As can be seen from the figures, neither mutual intersection nor gap between the yarns is 



Experimental methods for evaluating and quantifying geometric features generally fall into two clas- 
sifications: destructive and non-destructive. Work has been carried out in both of these activities, 
and a brief description of the techniques and typical results follows. 


The specimen preparation is carried out by mounting a small sample (1.0 x 1.0 x sample thickness) 
in a matrix cylinder (2.0 cm in diameter and 2.0 cm in height). Successive manual polishing 
operations are carried out to prepare the specimen for microscopical observation. The polishing 




pixels. Each pixel can be coded using the available gray scale range between black and white 






Modelling of the geometrical reinforcement requires some experimental activities to validate the 
work. The use of established conventional x-radiography is in principle attractive for this purpose 
but requires enhancement of the yarns since graphite and epoxy are both transparent to standard 
operating voltages. Investigations into two different methods show that enhancement can be ac- 
complished and yarn detection is possible. One approach attempts to image the entire yarn body 
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involved the doping of individual yarns with substances possessing high attenuation coefficients* 
This was aimed at imaging the entire yarn in width and length in a radiograph and in cross-section 
with tomography. The second approach examined the introduction of radiographically opaque 
tracer fibers into the graphite yarns. Although the actual yarn would not be imaged, the idea was 
to determine yarn direction on a global scale. 



Samples for this study consisted of six 12 K graphite yarns unidirectionally woven with a 
continuous T40 Kevlar fill. Epon 8132 resin with V-40 hardener was used for the matrix at a 1:2.5 
mix ratio. Individual yarns were treated before weaving. 
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were 80 kV and 4.5 mA for 45 second exposure. 




Experimental Neutron Tomography Results 
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Both the 40 and 37 AWG magnet wires in Samples 6 and 7 were easily distinguished. Though the 
yarns were not discerned, their paths and orientations were readily evident. Attractive features of 
the magnet wire included its ease of processing, size, ductility and low cost. The idea for using 



tracer fibers in composites has been investigated for some time though typically plies rather than 



One of the objectives of determining the geometric properties of textile reinforcements is to provide 
a means to understand the mechanical response of the composite formed from the fabric. Since the 
mechanical properties of a continuous fiber reinforced composite are clearly highly dependent upon 
geometric properties, it is necessary to understand the structure of the reinforcement in order to 
carry out any property predictions. 


Two basic approaches to this problem were employed. In one case the geometric data were 
used as input to a meshing algorithm creating a finite element model of the yarns and surrounding 
matrix. In the other the geometric data was reduced to simpler orientation parameters and a 
smearing method was employed. 
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FGM provides a quick method for constructing good predictions of elastic behavior. With additions 
to the approach, strain dependent tensile behavior can be modeled, such as that associated with 
metal matrix composites (plastic deformation of matrix [15]) and ceramic matrix composites (micro- 
cracking of matrix [16]). 



The basic idea behind the FGM is to treat the fibers and matrix as a set of composite rods 



where C g i 0 bal and Ci OC ai are the global and local stiffness matrices respectively and T a and 
T e are the stress and strain transformation matrices successively. 
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The formation of the transformation matrices depends exclusively on the direction cosines 
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Through the use of geometric relationships, it is possible to determine the other two vectors. 
Since rl, r“2, and f 3 are mutually orthogonal, 
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Since l\ is the direction cosine of the fiber with the x-axis (cos(l,x)), the definition of mi 
corresponds with mi = sin( 1, x). This places the 2-axis of the fiber in the plane made with the fiber 
and the x-axis of the composite. The remaining 5 terms can be solved from 8 and 9. Of course, 
it is important to use the relationships in 8 for the determination of the last term in each vector 
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Then, for each system of reinforcing fibers, the global stiffness C g i 0 b a ii is calculated and the 
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